Theory of the metamagnetic crossover in CeRu2Si2 
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Based on the periodic Anderson model, it is shown that the competition between the quenching 
of magnetic moments by local quantum spin fluctuations and a magnetic exchange interaction caused 
by the virtual exchange of pair excitations of quasiparticles in spin channels is responsible for the 
metamagnetic crossover in CeRu2Si2, cooperated with the electron- lattice interaction. The strength 
of the exchange interaction is proportional to the bandwidth of quasiparticles and its sign changes 
with increasing magnetizations; it is antiferromagnetic in the absence of magnetizations, whereas it 
is ferromagnetic in the metamagnetic crossover region. Experimental results of static quantities are 
well reproduced. 
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I. INTRODUCTION 

The metamagnetic crossover in CeRu2Si2 is an im- 
portant issue. The compound has a large electronic 
specific heat coefficient of 7 ~ 360 mJ/mol K^,EJ and 
it shows a sharp ,j^rease of magnetization at the field 



of Hm ^ 7.7 
magnetostrictiorcl'l 



Other physical, properties such as 
and specific heatu are also anomalous 
in this field region. One of the most crucial experimcr 
results to be explained is the single-parameter scaling;t 
independent experimental quantities for different pres- 
sures are scaled with a single energy parameter ksTx, 
where Tr- is called the Kondo temperature. The Kondo 
temperature is the energy scale of local quantum spin 
fluctuations and is approximately equal to the bandwidth 
of quasiparticles. P 

The present authors showed in a previous papem that 
an exchange interaction caused by the virtual exchange 
of pair excitations of quasiparticles in spin channels plays 
a critical role in the metamagnetic crossover as well as 
the electron-lattice interaction called the Kondo volume- 
collapse effecLu This exchange interaction has the fol- 
lowing two novel properties. First, it changes from being 
antiferromagnetic at zero fields to being ferromagnetic in 
the metamagnetic crossover region because of a pseudo- 
gap structure in the density of quasiparticle states, which 
is characteristic of heavy-electron compounds. Second, 
the strength of the exchange interaction is proportional 
to the bandwidth of quasiparticles. Then, the single- 
parameter scaling can be easily explained. 

After the previous paper was submitted, a detailed 
measurement of the field-dependence of specific heat was 
reported by Aoki et ali3 At sufficiently low tempera- 
tures, a single sharp peak was observed at Hm, while 
a double-peak structure was found at higher tempera- 
tures. They argued that such a result can be explained if 
a sharp peak exists in the density of states. A similar ar- 
gument applies straightforwardly to our previous model; 
however, if the pseudogap is enough deep, the magneti- 



zation process shows a first-order transition within the 
theoretical framework of the-jicevious paper. Recently, 
two theories were proposedO-Ej They claimed that an 
anisotropic c-f mixing plays an important role in the 
metamagnetic behavior. In both theories, however, the 
peak structure in the density of states is too sharp be- 
cause the k dependence of hybridization matrices is im- 
properly treated; one-dimensional van Hove singularity is 
irrelevant. Such an extremely sharp peak is inconsistent 
with the experimental result of specific heat. Further- 
more, it should give a discontinuous transition instead of 
the metamagnetic crossover if the novel exchange inter- 
action mentioned above is properly taken into account. 
To reproduce the experimental results of magnetization 
and specific heat simultaneously has not been achieved 
so far. 

A main purpose of this paper is to reformulate and 
improve the previous theory based on the periodic An- 
derson model so as to explain static properties of the 
system. We will also give a brief comment on intersite 
spin fluctuations around the zone center. 



II. MODEL 



The periodic Anderson model is written as 



^ = y]-E^A(k)4k^aAka 
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(2.1) 

with A the band index of conduction electrons, rii^ — 
liafiT, ^nd H* = rrioH, where mg is the saturation mag- 
netization per / electron. Other notations are standard. 
The kinetic energy of conduction electrons, ExCk), and 
the / electron level, Ef , are measured from the chemical 
potential, respectively. 



When multiple conduction bands are assumed, a pseu- 
dogap structure is expected in the density of states be- 
cause of the hvbridization between the / band and con- 
duction bands £3 Since the property of the magnetic ex- 
change interaction that plays a critical role in the mag- 
netization process depends on the shape of the density 
of quasiparticle states as shown in subsequent sections, a 
phenomenological model for the density of quasiparticle 
states is employed instead of considering explicit forms 
oi Ex(k) and VA(k) for simplicity. Furthermore, in order 
to make our treatment easy, we assume that the syst em 
is symmetrical. In other words, the Hamiltonian ( |2.l| ) is 
assumed to be invariant under the particle-hole transfor- 
mation. 

One of the most crucial issues in constructing a theory 
of strongly correlated electron systems is how properly lo- 
cal quantum spin fluctuations, which are responsible for 
the quenching of magnetic moments, are treated. They 
can be correctly taken into account in the single-site ap- 
proximation (SSA)llJ that is rigorous for paramagnetic 
states in infinite dimensions. Consider the zero-field case 
at first; H = 0. Within the SSA, Green's functions for / 
electrons and conduction electrons are respectively given 

by 



the MAM numerically. Therefore, Eq. (2.4) is a self- 
consistency condition for A(e). However, we do not per- 
form this self-consistent calculation in this paper. In- 
stead, several approximations will be taken in subsequent 
sections with the use of well-known results for the Kondo 
problem. 

Consider the spin susceptibility of the MAM, Xs{ii^i), 
where iuji is a boson energy with I an integer. The Kondo 
temperature for the periodic Anderson model, which is 
the energy scale of local quantum spin fluctuations, is 
defined by 



^imx.(+.0)^^-^. 



(2.6) 



In the same way as the previous paper, the electron- 
lattice interaction is taken into account simply through 
the volume dependence of the Kondo temperature: 



TKix)^TK{0)e~^, 



(2.7) 



with X = T Ay/Vb. Here, T - 190 is the Griineisen 
constant of Tk ■ For the sake of simplicity, the argument 
X will be omitted unless particularly required hereafter. 



G//^(ie„,k) 



iSn - Ef - ta{ien) - y^ 



|V3.(k)| 



iSn - Ex{k) 



(2.2) 



with Ect the single-site self-energy function, and 



GAA'cr(«en,k) = Sxygx{ien,'k-) 

+gxiiSn,'k)Vxik)Gff^iien,'k)V^,i'k)gx'{iSn,'k), (2.3) 

with gx{ien,k.) = [je„ — ExCk)]^^. Here, ie„ is an imagi- 
nary fermion energy with n an integer. To obtain the 
single-site self-energy function is, reduced to solving a 
single-impurity Anderson modelEj that has the same lo- 
calized electron leve l Ef and on-site Coulomb repulsion U 
as those in Eq. ( |2.l| ) . Call this Anderson model a mapped 
Anderson model (MAM) . Other parameters in the MAM 
are determined through the mapping condition: 



GffaiiEn) = — ^G//^(i£„,k), 

k 



(2.4) 



where N is the number of unit cells and Gffa is the 
Green's function of the MAM written as 



GffaiiSn 



1 



iSn - Ef - J^aiiSn) ~ L{ien) ' 



(2.5) 



with L(ie„) = (l/7r)/^de A(£)/(^£„ - e). Here, A(£) 
is the hybridization energy of the MAM. Once a trial 
function for A(£) is given, So.(i£n) is obtained by solving 



III. FERMI LIQUID DESCRIPTION 

The self-energy function is expanded as 

E,(z£„) = E,(+zO) + [1 - ,^,„]z£„ + ■ • • , (3.1) 

for small |£„|, where 0m is a mass enhancement factor in 
the SSA. Note that Eo-(-l-iO) = —Ef in the symmetrical 
case. Then the coherent part of Eq. (2.2) is written as 



.(c) 
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and correspondingly, G\l,^ is given by Eq. ( |2.3| ) with 

replacing Gffa by Gfi^. Quasiparticles are defined as 
the poles of Eq. (3.2), namely, the dispersion relation of 
quasiparticles is obtained by solving the following equa- 
tion 



4>mz - y^ 



I^A(k)P 

z-Ex{k) 



0. 



(3.3) 



We write the solutions as z = ^,y(k) with v representing 
the branch of q uasiparticles. By dividing the right-hand 
side of Eq. (3.2) into partial fractions, it follows that 



G*;]„(i£„,k) = y ^i 

^^"^ ' ^ Zl{k) ^e 



1 



e.(k)^ 



(3.4) 



where the renormalization factor is given by 



zl{]^) 



Y.'^x.i^) 



(3.5) 



with 7/A.(k) = |FA(k)|V[6(k) - ^A(k)]2. Similarly, we 
have 



Gil(*£n,k) = 5]-i 



1 



with 



Z,^(k) ze„-e.(k)^ 
?7A,.(k) 



^-^W ^™ + 5]77A.(k) 



(3.6) 



(3.7) 



It immediately follows from Eqs. (3.5) and (3.7) that the 
renormalization factors satisfy the relation 






(3.8) 

With the use of Eq. ( |3.8| ), the density of quasiparticle 
states can be written in the form 

(3.9) 

with e+ — e + iO. The energy scale for quasiparticles, 
Tq, is defined by 



p*{0) ^ 



1 



(3.10) 



4fcBTQ 
On the other hand, ^^(k) obeys the following condition 



^ zl(k) ^J 



(3.11) 



which can be proved by comparing Eqs. (3^) and (3.4) 
for the limiting case of |e„| — > +cx). It should be noted 
that Zl (k) ~ (j)m for the quasiparticle band that is closest 
to the Fermi level and ZlCk) ^ <j)m for other branches 
for a given k, that is, 

7 r 1 for \U^)\<2kBTQ 

^"' 2± < . (3.12) 

ZU^) [O for |^,(k)|»2fcBrQ 

Consider a spectral function of renormalized / electrons, 
p*r, defined by 






P}(^ 



(3.13) 



or equivalently 



<5(£-^.(k)). 



k zlik) 

Equation ( |3.12 ) tells that 

' p*ie) for \e\<2kBTQ 
for |e| > 2kBTQ 



(3.14) 



P}i^) 



(3.15) 



and it follows from Eq. ( |3.11 ) that p*f satisfies the nor- 
malization condition of /_ JJ£ p*f{s) = 1. 

The Luttinger's argumentil3 applies straightforwardly 
to the periodic Anderson model. Consider the number of 
electrons per unit cell: 

L^e+-""|G//,(ze„,k) + ^GAA.(*e„,k)|, 



NP 



(3.16) 



with P = l/ksT. At T = K, Eq. ( |3.16D can be trans- 
formed to 



de /9*(e), 



(3.17) 



which is an exact relation and is identical to the Lut- 
tinger's theorem. Let us derive an expression for the 
specific heat. Recently an exact expression fox. the en- 
tropy of an interacting system was establishedJlj For the 
periodic Anderson model, it is written as 



-E 



ttN 



de^M^lj2lmln[-g-,^{e^,k)] 



dT 



-Im In 



-G-fUe+,\.) 



Re Gffa (£+ , k) • Im E^r (eh 



(3.18) 



where /(e) = l/(e'3^ + 1). In Eq. ( |3l8| ), Gjf^ can 
be replaced by its coherent part because the derivative 
of /(e) is non-zero only for small |e| at low tempera- 
tures. In differentiating the entropy with respect to T, 
we introduce two assumptions; Im Yia{£+) = 0, and 
(j)m = T-independent, that is, quasiparticles have an in- 
finite lifetime and temperature-independent mass. Since 
the chemical potential does not change as a function of 
temperature in a symmetrical model, the specific heat at 
constant volume is given by 



Cv 



T 



dS_ 
df 



E 



where the identity Td'^f/dT'^ 
been made use of. 



dee^p*{e), (3.19) 

-d/de{edf/dT) has 



The zero-temperature limit of Eq. ( ^.ig| ) is given by 
Cv = iT with 7 = (27r2fc|/3)//(0). On the other 
hand, a combination of Eqs. (2^), ( |3.13 ) and ( |3.15 ) gives 
/9*(0) ~ p}(0) = ^„i/7rA(0). When A(e) is constant, one 
can prove thatQ 



1 



ttA AksTx ' 



(3.20) 



and then it foUows Tk — T q in t he SS A, where Tk and 
Tq are defined by Eqs.(|2.6|) and ( 3.10| ), respectively. We 
take this approximation in order to estimate the Kondo 
temperature. From the experimental value of 7, we have 
Tk{0) c^ 38 K. 

Let us consider the finite-field case; H ^ 0. In the 
previous paper Ja we constructed a perturbation method 
to derive a microscopic Landau free energy, one of 
whose independent variables is the magnetization m = 
'Ylia'^^Jiafi'y)- ^^ ^hc prcscncc of magnetizations, physi- 
cal quantities, such as self-energy and polarization func- 
tions, can be expressed as a function of m instead of the 
magnetic field H . The equation that defines quasiparti- 
cles becomes 



5Y^a[m) ~ y^ 



I^A(k)P 

z-^A(k) 



0, 



(3.21) 



with 5Y,a{m) a magnetic part of the self-energy. In 
Eq. (3.21), m depende*jce of the mass enhancement fac- 



tor has been ignored. tZI The solutions for Eq. (3.21) are 
denoted by ^i/cr(k, m). All arguments for the finite- field 
case can be developed in parallel with those for the zero- 
field case with replacing ^i, (k) by S,ua (k, ttj) . For example, 
(l>,n/Zl^{]i.,m) satisfies a similar property as Eq. (3.12). 
Therefore, it follows by comparing Eqs. (3.3) and (3.21) 
that 



^,,(k,m)~^,(k)-CTAS(TO), 



(3.22) 



wher e aA E(m) = — (5So-(?7i)/^m. With the use of 
Eq. ( 3.22 ), the Luttinger's theorem gives 



= ^0-/ de p*f{e + aAE{Tn)), 

„ J — 00 



(3.23) 



where the polarization of conduction electrons has been 
ignored. In Eq. (3.23), we have replaced p* with p^ 
because the right-hand side of Eq. ( 3.25 ) is a differ- 



ence between the contributions from up and down spin 
directions, only the low-energy part is relevant. From 
Eq. (3.23), the magnetic part of the self-energy, AE{m), 
can be determined as a function of m. 

A similar treatment for the MA M is also possible. Con- 
sider the coherent part of Eq. (2.5), which is defined 



by G^jj^iiSn) = i^/N)J2]^G^jj^{ien,'k-) and written as 
[(pmiSn — i'-'^H*^")]"^ ^t zero fields. In the 



//tr^*^" 



Gyl„iien,m) 



(i)^[ien + (jAEA{m)] - L('=)(ie„) 



(3.24) 



Note that AEA{m) is different from AE{m) because the 
MAM i s d etermined in the absence of magnetizations 
by Eq (2.4) and then magnetic one-body potentials are 
added to bpth the periodic Anderson model and MAM 
separatelyH In other words, Ai?^(m) is a single-site term 
even with respect to m, whereas AE{m) includes multi- 
site magnetizations. Because the magnetization in t 
MAM is also given by m to leading order in kBTx/Ul 
AEAifn) can be approximately evaluated from 



de 



-— Im 






(e+,m). (3.25) 



In case of a constant hybridization energy, Eq.(3.25) is a 
rigorous relation and is nothing but the Friedel sum rule. 
It should b e not ed tha t the evaluation of AE and AEa 
from Eqs. ( 3.23| ) and ( ^.25| ) is one of the most essential 
improvements. In the previous paper, we assumed that 
AE{'m) = AEA{m) = (4A:BTif/7r) tan(7rm/2). 

Before closing this section, it is helpful to mention the 
volume dependence of the physical quantities that have 
appeared above. First, it follows from Eqs. (p]^) and 
( |3.2Cl| ) that (j^m oc e^. Taking notice of Eq. (|t|), we have 
p*f{e]x) = e^pj(ee=^;0) and L(=)(ie„;a;) = L^(ie„e^; 0). 
Therefore, Eqs. ( |3.23D and ( |3.25D give AE{m) oc e""^ and 
AEAifn) oc e~^, respectively. 



IV. MAGNETIC EXCHANGE INTERACTION 

In this section, we study magnetic exchange interac- 
tions working between quasiparticles. Consider the mag- 
netic susceptibility to this end. In the presence of magne- 
tizations, spin and charge channels of the susceptibility 
are coupled with each other in general. In the symmet- 
rical case, however, they are separated and the magnetic 
susceptibility can be expressed as 



Xs(«w/,q,"i) 



2iis{iuji,ci,m) 
1 - C/7r3(iw/,q, to)' 



(4.1) 



with Tig — \ Yliaa' "'"''^o-cr' , whcrc TTo-cr' is the irreducible 
polarization function. Similarly, the magnetic suscepti- 
bility of the MAM is given by 



Xs(,JWi,TOJ 



27^s(^u;;,TO) 

1 - UTTs{iLUl,m) ' 



(4.2) 



with TTs the single-site polarization function for the spin 
channel. The static an^ zero-temperature limit of Xs is 
approximately given bj 



presence of magnetizations, it becomes 



I 



(4.3) 




FIG. 1. Diagram for -n^^i in the site representation. A 
thick line with an arrow stands for the Green's function. Over- 
counted contributions should be subtracted (see the text). 



Because Xs is also of order l/ksTK^ Eq. (H.lD can be 
rewritten in the form 



Xs(Jw;,q,m) = 



Xs(iw;,m) 



1 - 4/s(iu;;,q, m)xs(it^/,m)' 



(4.4) 



where /s(iwi,q, m) is the intersite exchange interaction, 
which is given by 

Isiiuji,q,m) == 2U^{Trs{iuJi,q,m) - 7rs(iw;, m)}, (4.5) 



to leading order in ksTx/U. Note that Eq. (4.4) is con- 
sistent with the physical picture of the Kondo lattice that 
local quantum spin fluctuations at each site are connected 
with one another by the intersite interaction. In order to 
calculate Ig, we consider a two-line diagram shown in 
Fig. P for the polarization part tTo-ct' , which remains fi- 
nite even in infinite dimensions for specific q's. In other 
words, the two- line diagram is a leading term with re- 
spect to 1/d for specific q's, with d the spatial dimen- 
sionality. In Fig. |l|j^ Xa-T is a three-point vertex function 
of the MAM and X^t denotes that of the periodic An- 
derson model of d — > -|-oo limit. Note that X^r is a 
single-site term, while A^r contains multi-site terms with 
respect to both U and m. It follows from Fig. |l]that 



C/'v-rr 



Is{iuji,q,m) = — o" X! ^^'^^ («£n , «£« 
X Xs{iSn,ien+iuJi,(i,m) 



iuJi,m) 



I — ^ Gffcr (ie„ , k, m)Gffc (ie„ + iu;i,'k + q, m) 

k 

- Gffa{ien,m)Gffa{ien + iuji,m)\, (4.6) 



where Xs = AL 



A|| and As 



^TT - ^iT' respec- 
tively. In Eq. (4.6), the second term is necessary not 
only to subtract the single-site portion but also to avoid 
overcountings. j— , 

The main contribution to Is is divided into two parts :t3 



Is {iuJi , q, m) = Js (q) -|- Jq {iui , q, m) . 



(4.7) 



The first term is an exchange interaction caused by 
the virtual exchange of high-energy spin excitations. 



which scarcely depends on m. It is composed of three 
terms; the conventional superexchange interaction, an 
extended superexchange interaction and an extended 
RKKY interactionlla Because the property of Js depends 
on the whole band structure, we treat it as a phenomeno- 
logical parameter. The second term, jQ{iLUi, q, m), is due 
to the virtual exchange of low-energy spin excitations 
within quasiparticle bands: 



jQiiuJi,q,m) 



El 



^A(^)(+iO,+iO,m) 



xXs{+iO,+iO,0,m) 

(4.8) 



■ Gy}^ (ie„ , m)Gyj^ (ie„ + iw/ , m) | , 



where energy dependence and q dependence of the vertex 
parts have been ignored. The vertex functions in Eq.( [l.§| ) 
are related to the magnetic parts of the corresponding 
self-energy functions: 



and 



U-UA), -r. -n N T dAEA{m) 
-A(■4)(+^O,+^O,m) = 0,„ g^^, 



U~ , ,7 dAE(m) 

-Xs{+iO, +iO, 0, m = (j)^ — -^-^ 
2 am 



(4.9) 



(4.10) 



respectively. Therefore we have 



jQii^iiq,^) 
_ OAEa dAE 
dm dm 

with 



[P{iuji,q,m)-P{i0ji,m)}, (4.11) 



P{iuji,q,m) 

^ J_ v^ /(^^'a(k + q,m))-/(^^^(k,TO)) 
N ^ ^^^(k, m) -^t.'o-(k-t-q, m) -f iw; 



(4.12) 



where Eq. ( |3.12D has been used and 

P(iw;,TO) 



^ 
/? 



^ &f%{ien, m)G)%{ie„ + iLUi,m). (4.13) 



In the following part of this se ction^ we study the static 
and uniform component of Eq. ( 4.1l| ), Jq (-|-iO, 0,m). It 
is abbreviated to Jq^m, x) with the argument x explicitly 
written. Because excitations between different branches 



we obtain 



in Eq. (4.12) can be neglected for low-energy phenomena, 
in 

P{+iO, 0, m) = Y^ p}{aAE{m)). (4.14) 



p/(e) 



c, =0.17 
Ci = 0.22 





FIG. 2. Phenomenological model for P/(e) for three values 
of ci. Energies are in units of ksTK- 



On the other hand, Eq. ( |4.13| ) is expanded as 



P{uj+,m) 



dAEA 



dm 



TT 
— ILU 

2 



OAEa 
dm 



(4.15) 



with the use of Eq. (1 3.25^, w he re uj^l - — uj -MO. It can be 
easily shown from Eqs. ( |4.1l| ), ( [1.14 ) and ( |4.15| ) that the 
volume dependence of Jq (m, x) is the same as that of Tpc , 
namely, jQ{m,x) — e^^ jQ{m,Q) and the magnitude of 
Jqim,, x) is of order UbTk- Thus, Jgim,, x) is scaled with 
Tk, or the bandwidth of quasiparticles. This_property is 
responsible for the single-parameter scaling.olj Once an 
explicit form of p*f{s) is given, jQ(m, 0) is straightfor- 
wardly calculated as a function of to. As a phenomeno- 
logical model for p*f{s), we employ the following model 
in the same way as the previous paper. 



P/(£) = -r^^{-K-[Die;c2,C3)+D{e;-C2,C3)] 

kBJ-K ^ 2 



+ (l-ci)i?(e-;0,C4)}, 



(4.16) 



with nD{y; a, b) — b/[{y — a)^ + b^] and e = e/kBTx- In 
this paper, C2 = 0.60 and C3 = 0.15 are assumed, which 
are different from previous ones. Because of the condi- 
tion p*f{0) ~ I/AUbTk-, only ci is variable. Figure || rep- 
resents p*f{e) for three values of ci, and the corresponding 
results for Jq(to,,0) are shown in Fig. y. The exchange 
interaction Jq changes from being antiferromagnetic to 
being ferromagnetic with increasing m, which is consis- 
tent with the previous result. The maximum value of Jq 
becomes larger when we raise the height of the peaks in 
p*f. Furthermore, it can be checked by changing C2 with 
Ci fixed that shifting of the location of the peaks to the 
band-edge side also enhances ferromagnetic instability. 



FIG. 3. Exchange interaction Jqira, 0) (in units of ksTK) 
as a function of m. 



V. THERMODYNAMIC QUANTITIES 

The thermodynamic potential is expressed as 

fl{m, x) = f^para(a^) + ^(to,, x) + AJ7(rn, x). (5.1) 

Here, the first term is that for the paramagnetic state 
including a contribution from the lattice system. It is 
written asl 



^para(a;) = fipara(O) + fcsTXO) 



2k, 



(5.2) 



where k is a dimensionless compressibility. The second 
term and the last term in Q are the magnetic single-site 
and multisite terms; they are defined by d^Cl/dm^ = 
l/Xs{+iO,m), and d'^ATL/dm^ = -Is{+iO,0,m)/A, re- 
spectively. By integrating them, we obtain 



^{m,,x) = kBTK{x){\ — \J\ — m? ), 



(5.3) 



and 



AQ.{m, x) 

Js(0) 

— ) 

8 



dm 
Jo 



dm ' 



jQ{m" ,x) 



(5.4) 



In this section, we confine our study to the case of ci = 
0.17. The equilibrium values of to and x are determined 
by solving the simultaneous equations, dQ/dm = H* and 
dQ/dx — 0; they are shown in Fig. ^. Here, we have set 
the values of three parameters as Js(0) — 0.3fc_BTft:(0), 
K — 2.4, and mo = 2.0/iB m order to have a better fit 
with experiments. The value of k corresponds to 1.09 
[Mbar]~^ in the ordinary unit. It should be mentioned 
that if the magnetization process is calculated with vol- 
ume fixed constant against the magnetic field, J.he meta- 
magnetic crossover is considerably suppressed. B 




H [T 



FIG. 4. Magnetization and magnetostriction. Ejajwrimen- 
tal data are shown by dots and circles, respectively.crEl 



Let us next study the specific heat. In the pres- 
ence of magnetizations, it is given by Eq.( |3.19[ ) with 
p*{e) replaced by p*A£ + alS.E{m)). Figure g shows 
Cv /T — H curves for three values of temperature; T = 
K, 1.5 K, and 3 K. The double-peak structure is quali- 
tatively reucpduced in the present model. Note that the 
experimentllj was performed at constant pressure. From 
the experimental data, the. specific heat at constant vol- 
ume has been estimated;t2l it is also shown in Fig. || for 
T = 1.5 K. In our calculation, no double-peak appears at 
T — 1.5 K, whereas it is observed above T ~ 1.0 K in the 
experiment. A probable reason for this discrepancy is th e 
easy estimation of the value of Tk based on Eq. (3.20). 
In addition to that, a contribution from ferromagnetic 
spin fluctuations or paramagnons, which is excluded in 
the SSA, may not be negligible. Figure ^ shows the tem- 
perature dependence of Cy jT for various values of the 
field, which is consistent with the experiment (Fig. 2 of 
reference 6). 



FIG. 5. Cv jT as a function of H . Experimentsidata at 
constant pressure for T = 0.25K are shown by dotsJlj Circles 
represent the valuffl&iof Cv jT at T = 1.5 K evaluated from 
experimental data.tJ 
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FIG. 6. Cv It versus T for various values of the magnetic 
field. 



Finally, the isothermal compressibility is given by 



1 



t^riii) 



^0 






dxdm 



dm 
dx 



H 



(5.5) 



As shown ip-Fig. R, it is strongly enhanced at Hm- Mat- 
suhira et atB have studied the field variation of the com- 
pressibility by using thermodynamic relations and assum- 
ing the single-parameter scaling. Our calculation shows 
a good agreement with their result. 



VI. DISCUSSION 
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It is worth while to comment on spin fluctuations 
around the zone center, q ~ 0. The imaginary part of 



Eq. (4.12) is expanded as 



Im P(LUjf., q, m) ~ c{m) ■ h 
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(6.1) 



H [T 
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FIG. 7. Isothermal compressibility as a function of H. 

for small q and small ti^/|q|, where the coefficient c{'m) 
depends on the band structure. On the other hMd, the 
inverse of the local susceptibility is expanded asc3 



Xs{^ 



'+>' 



Xs(+iO, m) 



-lUJ - 



(6.2) 



Because the w-linear terms in Eqs. (4.15) and ( |6.2|) can 
be neglected compared with the uj/q term in Eq. ( |6.lD for 
small q,EJ the total susceptibility is written in the form 

X6(w+,q,m) 



Re [l/xs(w+,q,TO)] -iC{m) -uj/q 

(6.3) 

with C{ra) = d/S.EA/ dm-d/S.E /dra-c{ra). Since the sys- 
tem is close to ferromagnetic instability in the metamag- 
netic crossover region; Re_[l/xs(+«0, 0, m)] <C ksTj^, the 
imaginary part of Eq. ( |6.3| ) is expected to be enhanced 
around the metamagnetic point for small uj. 

On the other hand, spin fluctuations in the Kondo 
lattice have been studied assuming the w-independent 
RKKY interaction, 



Xs(w+,q,m) 



Xs{^+,m) 



'RKKY 



{n)XsiuJ+.m)'' 



(6.4) 



instead of Eq. (4.4) without any justification. In the low- 
energy region, however, Eq. (6.4) is expressed as 



Xs(w+,q,m) 



1 



Re [l/xs{i^+,(l,m)] -TTiuj/2' 



(6.5) 



and behaves quite differently from Eq. ( |6.3[ ) arou nd the 
zone center. Thus the argument based on Eq. (6.4) is 
apparently inadequate. Indeed, recent results of neutron 
scatt erin g measurementscj cannot be understood from 
Eq. (F 
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